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Abstract A measurement of the top quark pair pro- 
duction cross section in the final state with a hadron- 
ically decaying tau lepton and jets is presented. The 
analysis is based on proton-proton collision data re- 
corded by the ATLAS experiment at the LHC, with a 
centre-of-mass energy of 7 TeV. The data sample cor- 
responds to an integrated luminosity of 1.67 fb _1 . The 
cross section is measured to be o t i — 194 ± 18 (stat.) ± 
46 (syst.) pb and is in agreement with other measure- 
ments and with the Standard Model prediction. 

1 Introduction 

Top quark pairs (tt) are produced in abundance at the 
Large Hadron Collider (LHC) due to the high centre- 
of-mass energy of 7 TeV. The large sample of tt events 
collected with the ATLAS detector makes it possible to 
study experimentally challenging decay channels and 
topologies. This letter describes a measurement of the 
tt production cross section. The final state studied here 
consists of a hadronically decaying tau lepton (rh a d) 
and jets, corresponding to the tt — > [^had^T ][&<?<?] decay, 
where b and q are used to denote 6-quarks and lighter 
quarks, respectively. Such an event topology with a 
hadronically decaying tau lepton corresponds to ap- 
proximately 10% of all tt decays [l|. 

A tt cross-section measurement in the final state 
with tau leptons makes it possible to probe flavour- 
dependent effects in top quark decays. It is also rel- 
evant to searches for processes beyond the Standard 
Model, where tt events with tau leptons in the final 
state are a dominant background. This measurement is 
particularly important for hypothetical charged Higgs 
boson production hi top quark decays, where the 
existence of a charged Higgs boson would lead to an 



enhancement in the cross section for the considered tt 
final state. The measurement presented here is com- 
plementary to the previously published tau + lepton 
(electron or muon) channel measurement Q . The most 
recent cross-section measurements of the tau + jets de- 
cay channel have been performed by the CDF and DO 
collaborations in proton-antiproton collisions at i/s = 
1.96 TeV l7j, 8]. This is the first measurement reported 
in this specific channel at the LHC. 

In this analysis, events with at least five jets are 
selected, where two of the jets are identified as hav- 
ing originated from &-quarks. After identifying the two 
jets likely to come from the hadronic decay of one of 
the top quarks, one of the remaining jets is selected 
as the ihad candidate from the other top quark. The 
Thud contribution is separated from quark- or gluon- 
initiated jets with a one-dimensional fit to the distribu- 
tion of the number of tracks associated with the Th a d 
candidate. Since the T^ad decays preferentially to one or 
three charged particles (and other neutral decay prod- 
ucts), this variable provides good separation between 
hadronically decaying tau leptons and jets, as the latter 
typically produce a large number of charged particles. 
The main backgrounds to the tt signal are multijet ev- 
ents, tt events with a different final state or signal events 
where the wrong jet is chosen as the Th a d candidate. A 
small contribution from single-top and W+jets events 
is also present. The distributions for the backgrounds 
used in the fit are obtained with data-driven methods. 



2 The ATLAS detector 

The ATLAS detector [9j is a multipurpose particle phy- 
sics detector with a forward-backward symmetric cylin- 
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drical geometry and a near-47r coverage in solid ang 
The inner tracking system covers the pseudorapidity 
range 1 77 1 < 2.5, and consists of a silicon pixel detec- 
tor, a silicon microstrip detector (SCT), and, for 1 77 1 
< 2.0, a transition radiation tracker. The inner detec- 
tor is surrounded by a thin superconducting solenoid 
providing a 2 T magnetic field along the beam direc- 
tion. A high-granularity liquid-argon sampling electro- 
magnetic calorimeter covers the region 1 77 1 < 3.2. An 
iron/scintillator tile hadronic calorimeter provides cov- 
erage in the range | 77 \ < 1.7. The end-cap and forward 
regions, spanning 1.5 < | rj \ < 4.9, are instrumented 
with liquid-argon calorimeters for both electromagnetic 
and hadronic measurements. The muon spectrometer 
surrounds the calorimeters. It consists of three large 
air-core superconducting toroid systems and separate 
trigger and high-precision tracking chambers providing 
accurate muon tracking for 1 7/ 1 < 2.7. 



and as uncertainties according to the MSTW prescrip- 
tion [15j, and cross-checked with the next-to-leading- 
order -I- next-to-next-to-leading-log calculation of Cac- 
ciari et al. [l6| as implemented in TOP++ 1.0 17|. Tau 
lepton decays are modelled with Tauola ■ Samples 
of simulated events are also used to estimate the small 
contributions from H^+jets, Z+jets, single-top-quark 
and diboson events, as described in Ref. [l9(. The gen- 
erated events were processed through the full ATLAS 
detector simulation using Geant4 2^, 21 1, followed by 
the trigger and offline reconstruction. The distribution 
of the number of pile-up events (i.e. collisions in the 
same, or nearby, bunch crossing as the hard-scattering 
event) is adjusted to match the scattering multiplicity 
measured in the data. 



4 Event selection 



3 Data and simulation samples 

The data used in this analysis were collected during the 
first half of the 2011 data-taking period and correspond 
to a total integrated luminosity of £ = 1.67 fb _1 . The 
data sample was selected with a 6-jet trigger that re- 
quired at least four jets identified with | 77 | < 3.2 and a 
transverse energy (Et) above 10 GeV. Two of these jets 
were required to be identified as 6-jets using a dedicated 
high-level-trigger 6-tagging algorithm [l^] • This trigger 
was enabled for only part of the 2011 data-taking pe- 
riod and is therefore the limiting factor in determining 
the integrated luminosity of the dataset used. 

The selection efficiency for the ti —> Thad + jets 
signal is derived from Monte Carlo (MC) simulations. 
The MC@NLO v4.01 [HI generator, with the parton 
distribution function (PDF) set CT10 [l3|, is used for 
the ti signal. The theoretical prediction of the ti cross 
section for proton-proton collisions at a centre-of-mass 
energy of yfs = 7 TeV is a t i = 1672^ pb for a top quark 
mass of 172.5 GeV. It has been calculated at approx- 
imate next-to-next-to-leading order (NNLO) in QCD 
with Hathor 1.2 [HI using the MSTW2008 90% con- 
fidence level NNLO PDF sets [11], incorporating PDF 



1 ATLAS uses a right-handed coordinate system with its ori- 
gin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The z-axis 
points from the IP to the centre of the LHC ring, and the 
j/-axis points upward. Cylindrical coordinates (r, (f>) are used 
in the transverse (x, y) plane, tj> being the azimuthal angle 
around the beam pipe. The pseudorapidity is defined in terms 
of the polar angle 8 as 77 = — lntan(0/2). The variable AR is 
used to evaluate the distance between objects, and is defined 
as AR = <y(A<t>) 2 + [Ar]) 2 . 



Jets are reconstructed from clusters of calorimeter 
cells (iif using the anti-fc 4 algorithm 23, 24 1 with a 
radius parameter R = 0.4. The jets are calibrated us- 
ing transverse momentum- and 77-dependent corrections 
obtained from simulation and validated with collision 
data [25| . Candidate events are required to contain at 
least five jets with a transverse momentum (pt) larger 
than 20 GeV and 1 77 1 < 2.5. 

The identification of jets originating from fe-quarks 
is performed using algorithms that combine secondary 
vertex properties and track impact parameters 26] . The 
algorithm identifies 6-jets with an average efficiency of 
60% and provides a light-quark jet rejection factor of 
about 340 in ti topologies. The likelihood of misidenti- 
fying a Thad as a 6-jet in a ti event is approximately 5%. 
The two jets with the highest 6-tag probability are cho- 
sen as the event ^-candidates; events with fewer than 
two 6-jets are rejected. 

The magnitude of the missing transverse momen- 
tum (E 1 ™ 1SS ) is reconstructed from energy clusters in the 
calorimeters. The calibration of each cluster depends 
on the type of physical object associated with the clus- 
ter. The transverse momentum of muons in the event is 
also taken into account. The E^ 183 significance (S E m.i„ e ) 



is defined as E™ iss /(0.5 [%/GeV] • y/EZh), where £ E T 
is the scalar sum of the transverse momentum of all 
objects. Using a S E mis a requirement instead of a direct 
E™ lss requirement allows the rejection of multijet events 
where the £"™ ISS arises from energy resolution effects, 
while still retaining high efficiency for signal events with 
E™ lss coming from particles which do not interact with 
the detector [27J ■ Candidate events are required to have 

S E mias > 8. 
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Events containing a reconstructed electron or 
muon [H, [H with p T > 15 GeV and 1 r\ \ < 2.5 are 
vetoed to reduce the background due to events contain- 
ing W bosons that decay to electrons or muons and to 
avoid overlap with other ti cross-section measurements. 

In each event, a single Thad candidate is selected 
from the reconstructed jets using the following proce- 
dure. First, the reconstruction of the hadronically de- 
caying top quark is attempted by selecting the three jets 
(including exactly one of the two 6-candidates) which, 
when their four-momenta are added vectorially, give the 
highest px sum. The remaining jet with the highest px, 
excluding the remaining ^-candidate, is selected as the 
Thad candidate. Events where the Thad candidate px is 
below 40 GeV are rejected. 

The main contributions to the selected Thad candi- 
dates in the signal region come from the signal (rh a d 
from events), electrons from ti events and misidenti- 
fied jets from tt, single-top-quark production, VF+jets 
and multijet events. The contributions from Z / "/* +jets 
and diboson processes are negligible. 

5 Data analysis 

The majority of ri ia d decays are characterised by the 
presence of one or three charged hadrons in the fi- 
nal state, which can be reconstructed as charged parti- 
cle tracks in the inner detector. The number of tracks 
(?Hrack) originating from the interaction point associ- 
ated with a Thad candidate is used to separate the Thad 
contribution from the misidentified jet background. 
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Fig. 1 Distribution of n tra ck for Thad from MC ti events 
(solid black line), electrons from MC tt events (dashed red 
line), and for jets from multijet events from data (blue tri- 
angles). The multijet event selection uses a S^miss sideband 
region as described in Sect. [5] All distributions are normalised 
to unity. 



All selected tracks with px > 1 GeV located in a 
core region spanning AR < 0.2 around the jet axis are 




counted. To increase the discriminating power, tracks in 
the outer cone 0.2 < AR < 0.6 are also counted, using a 
variable px requirement that is dependent on both the 
AR of the outer track and the pt of the core tracks. 
This variable px requirement is designed to reduce the 
contribution from pile-up and underlying event tracks, 
and is explained in Ref. [3o| . The separation power of 
the ntrack variable is illustrated in Fig. [T] where a com- 
parison of the ntrack distribution is shown for Thad) elec- 
trons and misidentified jets from multijet events. 

To extract the signal from the n tra ck distribution, 
the data sample is fitted with three probability density 
functions (templates): a tau/ 'electron template, a gluon- 
jet template and a quark-jet template. The Thad com- 
ponent from tt events constitutes the signal in the event 
sample. Real electrons from tt events (either prompt or 
from leptonic tau decays) which failed to be rejected 
by the electron veto also contribute significantly to the 
event sample. The electron and Thad templates are com- 
bined into a single tau/ electron template to ensure a 
stable fit, using MC predictions to determine their rel- 
ative contributions. The tau/ electron template is ob- 
tained from simulated tt events. The small expected 
contributions to the real tau/electron component of the 
fit from single-top-quark and W+jets events do not 
change the shape of the template. 

The remaining significant contributions come from 
misidentified jets, and are separated into two templates. 
The gluon-jet template describes the QCD multijet pro- 
cesses which are dominated by gluon-initiated jets, and 
the quark-jet template describes the remaining processes 
(ti, single-top quark and VF+jets) that are enriched in 
quark-initiated jets. 

The gluon-jet template is determined using a side- 
band region where the S E ^ia S requirement is changed 
to 3 < S E! niss < 4. This selection greatly enhances the 
contribution from multijet events, reducing other con- 
tributions (e.g. from ti events) to less than 1%. The 
regions defined by the selection 2 < S^mi^ < 3 and 4 



< S f 



< 5 are also used to study any correlations 



between the S^miss criteria and the ntrack distribution. 

The quark-jet template is obtained from a ti control 
sample where the Thad candidate is replaced by a muon 
candidate. The reconstructed muon 29] is required to 
have px > 20 GeV, 1 77 1 < 2.5 and no jet within a dis- 
tance AR = 0.4. The requirement on the number of 
non-6-tagged jets is changed from three to two as the jet 
corresponding to the Thad is now replaced by a muon. 
The other selection requirements are the same as for 
the signal region. This isolates ti events with very high 
purity; the contribution from backgrounds is estimated 
from MC predictions to be at the 5% level, and consists 
mainly of single-top-quark and VK+jets events. The two 
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highest-px jets that are not identified as 6-jet candi- 
dates are selected as T ha d candidates. The template is 
corrected using MC simulations for differences in the 
transverse momentum distribution between the signal 
region and the control sample, and for the expected con- 
tribution to the control sample from tt dilepton events 
(tt -> n + T had + X, tt -> \i + e + X). 

6 Results 

An extended binned-likelihood fit is used to extract the 
different contributions from the n tra ck distribution. To 
improve the fit stability, a soft constraint is applied 
to the ratio of quark-jet events to tau/ 'electron events, 
which are dominated by the same process (tt events). 
The constraint, based on MC predictions, is a Gaussian 
with a width of 19% of its central value. This width was 
estimated based on studies of the associated systematic 
uncertainties using the same methodology as described 
in Sect. [71 The statistical uncertainties on the fit param- 
eters are calculated using the shape of the fit likelihood. 
The systematic uncertainties on the shapes of the tem- 
plates are propagated using a pseudo-experiment ap- 
proach, taking into account the bin-by-bin correlations. 
This yields a final number of tau/ electron events of 270 
± 24 (stat.) ± 11 (syst.). 




Fig. 2 The n trac k distribution for Thad candidates after all 
selection cuts. The black points correspond to data, while 
the solid black line is the result of the fit. The red (dashed), 
blue (dotted) and magenta (dash-dotted) histograms show 
the fitted contributions from the tau/electron signal, and the 
gluon-jet and quark-jet backgrounds, respectively. 

The fit results are shown in Fig. [5J A comparison 
between the fit results, and the expected event yields 
from the MC predictions is presented in Table [T] The 
numbers are in good agreement. 

To extract the number of signal events, predictions 
from simulation are used to subtract the backgrounds 



Source 


Number of events 


tau/electron 


tt (Thad) 


170 ± 40 


tt (electrons) 


47 ± 11 


Single top 


12 ± 2 


W+jets 


9 ± 5 


Total expected 


240 ± 50 


Fit result 


270 ± 24 (stat.) ± 11 (syst.) 


quark-jet 


tt (jets) 


540 ± 160 


Single top 


24 ± 4 


W+jets 


21 ± 12 


Total expected 


580 ± 160 


Fit result 


520 ± 97 (stat.) ± 78 (syst.) 


gluon-jet 


Fit result 


960 ± 77 (stat.) ± 74 (syst.) 



Table 1 Comparison of the numbers of events from MC ex- 
pectations and from the results of the fit to the data for the 
three templates. The uncertainties on the MC expectations 
include the systematic uncertainties of the selection efficiency 
described in Sect. [7] No MC predictions are available for the 
gluon-jet contribution. 



from W+jets and single-top events (9 ± 5 and 12 ± 2, 
respectively) from the fitted number of tau/electron ev- 
ents. The number is then scaled by the expected ratio, 
N T /(N T + N e ), of Thad and electrons passing the selec- 
tion in the tt sample. This ratio is estimated from MC 
simulation to be 0.78 ± 0.03 (stat.) ± 0.03 (syst.). This 
yields a final number of observed signal events of N T = 
194 ± 18 (stat.) ± 11 (syst.). 

The cross section is obtained using a t t — N T / (C ■ e). 
The efficiency (e) is estimated from MC simulation to 
be (6.0 ± 1.4) xl0~ 4 . It includes the branching frac- 
tions for the various tt decays and the acceptance, and 
assumes Br(tt -^T had + jets) to be 0.098 ± 0.002 [l|. 
The efficiency is corrected for a 13% difference between 
MC simulation and data in the trigger and ^-tagging 
efficiencies [26| . The method used for obtaining the un- 
certainty on the cross section is detailed in Sect. [7J 

The cross section is measured to be a t i = 194 ± 
18 (stat.) ± 46 (syst.) pb. 



7 Systematic uncertainties 

A summary of all systematic uncertainties on the cross 
section is given in Table [5] 

The uncertainty on the selection efficiency due to 
the choice of the configuration for the MC simulation is 
estimated by using alternative MC samples and reweight- 
ing procedures. The difference in the efficiency obtained 
from various configurations is taken as the uncertainty. 
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Source of uncertainty 


Relative uncertainty 


ISR/FSR 


15% 


Event generator 


11% 


Hadronisation model 


6% 


PDFs 


2% 


Pile-up 


1% 


6-jet tagging efficiency 


9% 


Jet energy scale 


5% 


_E™ lsa significance mismodelling 


5% 


6-jet trigger efficiency 


3% 


Jet energy resolution 


2% 



Fit systematic uncertainties 



4% 



Luminosity 



4% 



Total uncertainty 



24% 



Table 2 Systematic uncertainties on the tt cross section. 



The uncertainty on the modeling of the ISR/FSR is 
taken into account by using AcerMC [3l| samples with 
specific tunes aimed at conservatively varying the 
amount of parton showering [jjj ]. The uncertainty due 
to the choice of the matrix element event generator 
is estimated by comparing tt samples generated using 
MC@NLO, Powheg dlHil, and Alpgen 0. To 
study the impact of different hadronization models, ev- 
ents generated using Powheg are processed with two 
different hadronization programs: Herwig and Pythia. 
The uncertainty of the choice of PDFs is estimated us- 
ing a number of current PDF sets 37 1. 

Uncertainties on the simulation of the detector re- 
sponse are taken into account using dedicated studies 
of the reconstructed physics objects (electrons, muons, 
jets, _E™ 1SS ). The uncertainties considered are associ- 
ated with the jet energy scale, jet energy resolution, 
^-tagging efficiency, trigger efficiency and the i?™ ISS cal- 
culation 2^, 26} . The uncertainty due to mismodelling 
of the lepton veto is estimated using the uncertainties 
on the muon and electron reconstruction efficiencies de- 
termined from independent data samples, and is found 
to be negligible. 

To obtain the uncertainty on the fit results, varia- 
tions are applied to the templates to describe various 
systematic effects. As the tau / 'electron template is taken 
directly from MC-simulated tt events, the systematic 
uncertainties on this template are taken from estimates 
of the mismodelling of the simulation. The dominant 
contributions come from variations in the amount of 
ISR/FSR in the simulation (1%), the modelling of the 
pile-up (I %), and the statistical uncertainties (1%). Un- 
certainties on the track reconstruction efficiency, jet en- 
ergy scale, and the ratio of Thad to electrons are found to 
be negligible. The quark-jet template is obtained from a 
fi + jets control sample of tt events in data. The dom- 



inant contributions to the uncertainty come from the 
statistical uncertainties (4%), the difference in shape 
between the [i + jets template and the expected quark- 
jet distribution, estimated from MC samples (2%), and 
the MC-based subtraction of the dilepton contribution 
(1%). The uncertainty on the MC-based kinematic cor- 
rection is found to be negligible. The gluon-jet template 
is derived from a background-dominated sideband re- 
gion with small values of S^mi^. The two sources of 
uncertainties are the dependence of the template on 
the ^miBs criterion of the control region, obtained by 
varying the S , E mi aB requirement (1%), and the statisti- 
cal uncertainty of the control region (1%). The total 
systematic uncertainty on the fit is found to be 4%. 

The uncertainty on the luminosity is calculated to 
be 3.7% as described in Ref. 38[. The total systematic 



uncertainty on the cross section is 24%. 



8 Conclusions 

This letter presents a measurement of the top quark 
pair production cross section in the final state corre- 
sponding to the tt — > [frThad^r] [bqq] decay. The mea- 
surement uses a dataset corresponding to an integrated 
luminosity of 1.67 fb _1 of proton-proton collision data 
at a centre-of-mass energy of 7 TeV recorded by the AT- 
LAS experiment at the LHC. The signal has been ex- 
tracted by fitting the number of tracks associated with 
tau lepton candidates using templates derived from sim- 
ulation for the tt signal and from the data for the back- 
grounds. 

The tt production cross section is measured to be 
a t t = 194 ± 18 (stat.) ± 46 (syst.) pb. This result is 
compatible with the highest precision ATLAS measure- 
ments 3l| 4^, and with the result of 186 ± 13 (stat.) 
± 20 (syst.) ± 7 (lum.) pb obtained in the complemen- 
tary tau + lepton (electron or muon) channel 



6|. It is 

also in good agreement with the theoretical prediction 
of 167±^ pb. 
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